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Classical bremsstrahlung of a massless scalar field $ is studied in gravity mediated ultra-relativistic 
collisions with impact parameter b of two massive point particles in the presence of d non-compact 
or toroidal extra dimensions. The spectral and angular distribution of the scalar radiation are 
analyzed, while the total emitted energy is found to be strongly enhanced by a d— dependent 
power of the Lorentz factor 7. The direct radiation amplitude from the accelerated particles is 
shown to interfere destructively (in the first two leading ultra-relativistic orders) with the one due 
to the $ — $ — graviton interaction in the frequency regime 7/6 < ui < j 2 /b in all dimensions. 

PACS numbers: 11.27.+d, 98.80.Cq, 98.80.-k, 95.30.Sf 



Contents 

1. Introduction and results 

2. The setup 

A. The action 

B. Equations of motion 

C. Iterative solution 

1. The formal expansion and zeroth order equations 

2. The first order equations 

3. The second order equation for $— radiation 

D. The solution for °$, h MN and l z M in M A x T d 

E. $— radiation - Basic formulae 

3. The radiation amplitude 

A. The local amplitude 

1. The 7 — > 00 limit. Mode and frequency cut-offs. 

B. The non-local amplitude 

C. The part of the radiation amplitude and destructive interference 

1. j™ in the frequency range u 3> 7/6 

2. j™ in the frequency range cj < 7/6 

D. The part of the amplitude 

E. Summary 

4. The emitted energy - Spectral and angular distribution 

A. The z— type component of radiation with oj ~ 7 2 /& 

B. The z 1 — type radiation with § ~ 1 

C. The cases d=0, 1 

D. The estimate of the zz' — interference part of radiation 

5. Summary of results 
Acknowledgements 

A. Notation 

1. KK mode decomposition and Fourier transformation - Notation and conventions 

2. Integration over angles and frequencies 

3. Useful kinematical formulae 

B. Destructive interference for 7/6 <uj < 7 2 /& 



2 



References 



1. INTRODUCTION AND RESULTS 

The idea of TeV scale gravity with large extra dimensions (LED) |l|-|4( has triggered a lot of activity in 
particle physics and gravitation theory. One of the most interesting predictions is the possibility of black hole 
production in colliders [5]. According to Thome's hoop conjecture (generalized to higher dimensions Q), for 
energies of colliding particles higher than the D-dimensional Planck mass M* (transplanckian regime) such 
black holes should be produced classically for impact parameters b < R$, where Rs is the Schwarzschild radius 
associated with the center-of-mass collision energy. To establish the creation of a black hole in the collision of 
ultrarelativistic particles one has to find closed trapped surfaces in the corresponding space-time. To this aim, 
an idea due to Penrose Q was put into the form of an elaborated model in D = 4 by D'Eath and Payne [§[ 
and with extra dimensions by Eardley and Giddings Q and further refined in [lol ] (see recent reviews [llj | and 
references therein). In that model one is after a solution of Einstein's equations with a special metric ansatz 
generalizing the Aichelburg-Sexl metric [l2j (see also [ljj]) . The ansatz amounts to replacing the gravitational 
field of two ultrarelativistic particles before the collision by colliding shock waves, while the collision region is 
described by some linear differential equation for a metric function, which is amenable to construct exact or 
approximate solutions. The closed trapped surface emerging in such a solution for appropriate initial energies 
and impact parameters leads to an estimate of the produced black hole mass and its difference from the initial 
energy is interpreted as the amount of gravitational radiation produced (see e.g. the recent paper [Tij and 
references therein). 

The colliding waves model (CWM) of black hole formation is certainly a very nice and perhaps the simplest 
possible one, designed to answer an intriguing question about the nature of transplanckian collisions. It gives the 
gravitational radiation loss for head-on and almost head-on collisions and also demonstrates that the black hole 
is indeed present in the collision region. However, we would like to discuss some subtle points concerning the 
applicability of this model to high energy particle scattering. In fact, various effects which were not taken into 
account in the simple version of CWM such as an extended nature of the colliding particles [l5j and their parton 
structure [16J, have already been discussed. It was shown by Meade and Randall [16| that taking the finite size 
into account leads to a substantial decrease of the cross-section predicted by the CWM. It was also emphasized 
that an even more critical effect is the radiative energy loss of the colliding partons before their energy is 
trapped inside a black hole horizon. The CWM seems to be able to clarify these issues, but in a closer look it 
probably cannot. In fact, the Aichelburg-Sexl solution is the limiting form of the linearized gravitational field 
of ultrarelativistic particle moving with a constant velocity. Therefore, presenting colliding particles as plane 
waves implicitly assumes that their radiation is negligibly small, otherwise the particle trajectories should be 
modified substantially by radiation friction (for a discussion of radiation reaction in extra dimensions see [l7jV 
Furthermore, the energy mismatch between the mass of the black hole and the initial energy of the colliding 
particles, interpreted in CWM as radiation loss, is found to be of the order of the initial energy Q, which is by 
no means small. 

Therefore, alternative methods of computing radiation losses in tr ansp lanckian collisions seem to be desirable. 
One such approach is the one by Amati, Ciafaloni and Veneziano [18[ relevant in four dimensions and based 
on the combination of string and quantum field theory techniques. Other recent work on this subject includes 
'nil using various analytical classical and semi-classical approaches. In the framework of purely classical D = 4 
general relativity, numerical simulations were also performed of collision of two scalar field balls interacting 
gravitationally via exact Einstein equations [20]. Gravitational radiation was extracted imposing appropriate 
boundary conditions. Gravitational radiation in collisions of higher-dimensional black holes (with non-compact 
extra dimensions) was studied numerically in (2lj . 

For the ultrarelativistic scattering in models with large extra dimensions a crucial question is whether radiation 
is enhanced due to the extended phase space associated with extra dimensions. It was argued by Mironov and 
Morozov [22j that in the case of synchrotron radiation the expected enhancement can be damped by beaming 
of radiation in the forward direction, suppressing the number of excited Kaluza-Klein modes. In the case of 
bremsstrahlung the situation is different, namely in [23| it was found, that the energy loss of ultrarelativistic 
particles under non-gravitational scattering at small angle contains an additional factor 7 d due to the emission 
of light massive KK modes. Qualitatively, this can be explained as follows. For non-gravitational scattering in 
flat space the impact parameter b, the radiation frequency lo and the angle of emission •d <C 1 with respect to 
the momentum of the fast particle in the rest frame of the other are related by 

ujb($ 2 +7 -2 ) < 1. (1.1) 

Frequencies near the cut-off frequency u> ~ 7 2 /6 are emitted in the narrow (D — 1) dimensional cone d < 7 _1 , 
while intermediate frequencies are emitted into a wider cone. The main contribution comes from the first region, 
and in this case the emitted momenta in directions transverse to the brane are of the order w/7. The number of 
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such light modes is of the order (Rcj/j) d ~ (R'y/b)' 1 (R being the size of the compact extra dimensions), giving 
an extra factor j d to the radiated power. Analogous enhancement was reported for gravitational bremsstrahlung 
in transplanckian collisions [23 ]. compatible with the numerical study of (2~l| . 

The purpose of this paper is to study scalar radiation in gravity-mediated collisions in the presence of extra 
dimensions. Our approach was initiated in 0, HI] where we have briefly formulated purely classical perturba- 
tion theory for transplanckian scattering in ADD valid for small scattering angles, or large impact parameters. 
Indeed, as emphasized long ago by many authors, most notably by 't Hooft [26| . the gravitational force not 
only becomes dominant at transplanckian energies, but gravity itself becomes classical, at least in some signif- 
icant impact parameter regime. Thus, the applicability of classical theory in the bremsstrahlung problem at 
transplanckian energies seems to have solid theoretical grounds. 

More specifically, our approach amounts to solving the two-body field-mediated problem iteratively. In elec- 
trodynamics this is a well known method, allowing to calculate spectral-angular distributions of bremsstrahlung 
in the classical range of frequencies small with respect to the particle energy haj <C E. In general relativity this 
approach was suggested in [27[ under the name of "fast-motion approximation scheme" and was further devel- 
oped and called "post-linear formalism" . It was applied to the gravitational bremsstrahlung most notably by 
Kovacs and Thorne (28|. We will use a momentum space version of this approach, developed in [2i|, which has 
the advantage that it allows a fully analytical treatment of the problem. In [23|,[24[ we extended this technique to 
models with extra dimensions either infinite or compact. In [231 ] we calculated scalar bremsstrahlung radiation 
in the collision of two ultrarelativistic point-like particles interacting via a scalar field in flat space-time. Here 
we consider the situation in which the particles interact gravitationally but emit scalar radiation. This is an 
intermediate step towards the full treatment of the gravitational bremsstrahlung presented briefly in [24| ■ The 
latter case has additional complications due to the tensor nature of the radiation field and it will be presented 
in full detail in a future publication. Here, we will assume that only one of the colliding particles is coupled 
to the scalar field, so that their interaction is purely gravitational. On the other hand, we will compute only 
the scalar radiation emitted by the system. The main novel feature is this case is that the system becomes 
non-linear due to scalar-scalar-graviton vertex. As a result, the effective source of radiation field in the flat 
space picture becomes non-local due to contribution of the field stresses. 

In four dimensions, it was shown long ago [30(, that the contribution from the high-frequency regime 7/6 < 
ui < 7 2 /& is suppressed by a {actor 7 -4 due to destructive interference between the local and the non-local 
amplitudes, the latter being due to the gravitational interaction of the mediating field. The remaining radiation 
is beamed inside the cone with angle I/7, it has characteristic frequencies u ~ 0(7/6) and emitted energy of 
order E ~ j 3 . There is in addition a sub- leading component of emitted radiation with frequencies u> ~ 0(j 2 /b), 
which is also beamed and has E ~ j 2 . 

In the higher dimensional case the situation is more complicated. The destructive interference is also present, 
but with growing d the relative contribution of u ~ jb increases faster, than that of ui ~ 7/6 due to competition 
of the angular integrals. 

The powers of 7 of the emitted radiation energy in all frequency and angular regimes in the presence of d 
extra dimensions are summarized in the Table below. 



d\ 


u> -C 7/6 


to ~ 7/6 


cj ~ 7 2 /& 


W > 7 2 /6 


I- 1 


negligible 
(phase space) 


£ d ~7 3 


E d ~ 7 rf+2 


negligible radiation 


1 


negligible 
(phase space) 


E d ~ 7 rf +! 


negligible radiation 


negligible radiation 



In the special cases of d = 1 or 2 extra dimensions there is an extra In 7 in the expression for the energy 
emitted in the regime (oj ~ 7/6, $ ~ 1/7)- The energy is measured in units of x^m' 2 f 2 /b 3d+3 , with ml the 
mass of the target particle, K£> the D— dimensional gravitational coupling, / the scalar coupling of m and b the 
impact parameter of the collision and the overall numerical coefficients can be found in the text. 
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2. THE SETUP 



A. The action 



The goal here is to calculate within the ADD scenario classical spin-zero bremsstrahlung in ultra-relativistic 
gravity- mediated scattering of two massive point particles m and ml. The space-time is assumed to be 
M4 x T d , the product of four-dimensional Minkowski space and a d-dimensional torus, with coordinates 
x M = (2^,2/), M = 0,1,..., D- l,n = 0,...,3, i = l,...d. 

Particles move in M4 (the brane) and interact via the gravitational field guN, which propagates in the whole 
space-time M4 x T d . We also assume the existence of a massless bulk scalar field $(x p ), which interacts with 
m, but not with ml . The action of the model is symbolically of the form 



S — Sn ~\~ S<£> ~\~ S m ~\~ S', 



m' ^ 



and explicitly, in an obvious correspondence, 



I? 



1 



,AIN 



■Al-N 

egMNZ z - 



(m+f®) 2 



dr 



lAl - IN 



e g M NZ z + 



dr' 

with 16ttGd = x?t) relating X£> to Newton's constant. Here the ein-beins e(r) and e'(r') are introduced, which 
lead to a somewhat unusual form of interaction with <£>, but it reduces to the standard non-derivative interaction 
once they are integrated out. The constant / is the scalar charge of m. To solve the corresponding coupled 
equations of motion we shall use perturbation theory with respect to the scalar and gravitational couplings. 
To zeroth order, both gravitational and scalar fields are absent and the two particles move along straight lines 
as determined by the initial conditions. In first order, one takes into account the (non-radiative Coulomb-like) 
gravitational and scalar fields produced by the particles in their zeroth order trajectories. Next, one computes 
the first order correction to their trajectories, due to their first order gravitational interaction (the scalar mutual 
interaction vanishes in our set-up). The leading contribution to the scalar radiation field, of interest here, emitted 
by the accelerated particle m is then obtained in the second order of perturbation theory. This approach allows 
us to compute consistently the lowest order radiation in the case of ultrarelativistic collision, when deviations 
from straight trajectories is small and the iterative solution of the coupled particle-field equations of motion is 
convergent. The resulting expression for the radiative energy loss will be therefore correct only in the leading 
ultrarelativistic order. 

To carry out such a computation it is sufficient to restrict oneself to linearized gravity. One writes gu n — 
Vmn + ^d^mn and replaces the Hilbert- Einstein action by its quadratic part 



Do 



-h MN a D h MN + -hU D h - -h MN d M d N h + -h MN d M d P h p N 



d D , 



(2.2) 



where the Minkowski metric is t]mn — diag(l, — 1, — 1, ...), Od = n 9m^n, raising/lowering the indices is 



performed with t]mn and h = hy^. To avoid classical renormalization (which in principle can be treated along 
the lines of [13]) we take into account only mutual gravitational interaction. At the linearized level the total 
gravitational field Iimn is the superposition of the fields h™ IN , h™ IN due to the particles m and m! , respectively, 
and of the gravitational field generated by the bulk scalar $. Assuming that the scalar interaction is of the 
same order as the gravitational one, the latter contribution to Hmn is °f higher order and will be neglected 
here. Therefore, to this order of approximation 



h-MN — h 



AIN 



l MN 



(2.3) 



For simplicity the superscripts m, m! will be omitted in what follows. Thus, h^jy and h\\ N will be denoted 
as Hmn an d h' MN , respectively. Ignoring the self-interaction and the associated radiation reaction problem, we 



will consider each of the particles ml and m as moving in the other's metric 

gMN = Vmn + x-ohuN and g' MN = rj M N + ><: D h' MN , 



(2.4) 



respectively. Correspondingly, the action for the particle m, which interacts with both gravity and takes the 
form 



Similarly, the action for ml is 



e(r) 



AIN 



AIN 



)i M. N+ (m+f^y 



e'(r,MN + XDhMN)z' M z' N + 



m ' 



dr. 



dr. 



(2.5) 



(2.6) 
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Finally, the action for the scalar field, which propagates in the gravitational field h' MN of the uncharged 
particle, expanded to linearized order of the gravitational field, is 

d D x. (2.7) 

In principle, $ propagates in the full gravitational field h mn 4- h' MN of both particles, but the singular product 
of Hmn and $ generated by the same point particle to corresponds again to the self-action problem, which is 
ignored here. The products of h' MN due to m! and $ due to to does not lead to singularities and correctly 
describe the situation. 



2 



hj r/MN — XDh' MN 



B. Equations of motion 



Varying the action with respect to z(t) and z'(r) one obtains the linearized geodesic equations of each mass 
moving in the gravitational field of the other: 

fa (eg MN z ) = -^9lr,m z z > fa V 9m nz ) = -^9lr,m z z ■ ( 2 - 8 ) 
Variation with respect to the einbeins gives 

-2 _ 9mn z z 1-2 9MNZ Z 

(m + /$) 2 m' 2 K 1 

Substituting this back into the particles' actions (|2.5H2.6[) one is led to their more familiar form 



S m = - I (m + /$)(< ? ; /7V i M i JV ) 1/2 rfr, S m , = -m! J (gMN^z'^dr, 

from which the scalar field equation is obtained 

M = -^fh'Dn^ + K D h' MN $< MN + fj {g l MN i M i N ) 1 'H D {x - z{t)) dr, (2.10) 

Note, once again, that only the gravitational field due to the uncharged particle m! enters this equation. 

Finally, the linearized Einstein equations for the metric deviation due to the two particles in the De Donder 
gauge 



dNh MN = 1 gM h 



2 

are obtained from (12.51 12.61) : 



WA _ _ / mMN MN T \ rr^AIN _ f„-M-N^ ( x ~ Z ( T )) 



D D h lvln = -K D ^T m " - r, m " jf^J , T MN = J 6 z z y ^ ^dr, (2.11) 
where T — and linearization of the metric factor is understood. Similarly, 

D D h' MN = -hd ( T' MN - V MN -^—) , T' MN = [ e'z' M z' N 6D ^^ T)) dr. (2.12) 



D-2J' J J—g 

To ensure that the particles move on the brane, it is enough to choose the initial conditions y l (0) — 0, y l {0) = 
and similarly for m! . Then, using the equations of motion, it is easy to check that the entire world- lines will 
be z M = z M (r), = z' M (r) and the energy-momentum tensors will only have brane components ,r''"' '. 
However, the metric deviations h MN and h will have in addition diagonal bulk components due to the trace 
terms in (f2TTj) and ([2~T2l . 



C. Iterative solution 



Even in linearized gravity the relativistic two-body problem can not be solved exactly, so one has to use some 
approximation scheme. With the particle masses to, to' taken of the same order and eventually equal, the model 
is characterized by three classical length parameters. Namely, the classical radius of the scalar charge [23j |: 




(2.13) 
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the D-dimensional gravitational radius of the mass m at rest 

r g = (x 2 D m)^ , (2.14) 
and the Schwarzschild radius of the black hole, associated with the collision energy y/s (24l |: 



1 

\/7T 



arm 

d+2 



(2.15) 



In the (initial) rest frame of the mass m! one has y/s = 2mm' j, where 7 = — v 2 is the Lorentz factor of 

the collision, v being the relative velocity of the colliding particles. So 

rs-tyy", "=2(dVl)- ( 2 - 16 ) 

It will be assumed that the parameters r g and Tf are of the same order, and both much smaller than the impact 
parameter b: 

r 9 ~r/<6 7 - 2!/ , (2.17) 

or, in terms of rs [24]: 

b-»r s Y- (2.18) 

Under this condition, as will be shown below, the deviation of the metric from unity in the rest frame of m! is 
small, i.e. xpfiMN z' M z' N <C 1, which is necessary for the validity of the present perturbative treatment. 

1. The formal expansion and zeroth order equations 
The next step is to solve these equations iteratively. For the particle-m world-line one writes 

Z M = Z M + i z M + _ ^ Z M = U M T + (2ig) 

where the order is denoted by a left superscript and to zeroth order the particle moves in Minkowski space-time 
with constant velocity u M , and with °z (0) = 6 another constant vector. Both vectors u , b M will be 
assumed to lie on the brane, i.e. to have only /i-components and, in addition, to be orthogonal b^u^—0. It will 
be shown that as a consequence of the equations of motion 1 z M also lies on the brane. However, it is convenient 
to keep Z?-dimensional notation in all intermediate steps. 
For the particle-m' one writes similarly 

z' M = °z' M + V M + . . . , °z' M = u' m t, (2.20) 

assuming that at r = the particle is at the origin. We choose to work in the rest frame of m', and specify the 
coordinate axes on the brane so that u'^ = (1,0,0,0), u M = 7(1,0,0, v), 7 = — v 2 , and & M = (0, &, 0,0), 

where b is the impact parameter. When needed, one may think of the brane-localized vectors as D-dimensional 
vectors with zero bulk components, e.g. u M = (u M , 0, . . . , 0). 
In a similar fashion, the bulk scalar is expanded formally as: 

$= °$+ !$ + .... (2.21) 

Substitute in (|2.10[) and set x£> = to obtain for the leading contribution to $ the equation 

D D °f = / J S D {x -ur-b) dr. (2.22) 

The ein-beins are also expanded 

e = ° e + \ + . . . , e ' = V + V + . . . . (2.23) 
According to (|2.9[) one obtains in zeroth order 

°e = m + f °$, °e' = m'. (2.24) 
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Finally, for the metrics one writes 

h>MN = °t l MN + h>MN + ■ ■ • , (2.25) 

and similarly for h' MN . The leading order contributions to the metrics are then obtained from Eqs. (|2.11[ [2~T2| 
with the zeroth order source on the right hand side, i.e. with 

°T MN =mj 5 D {x ~ °z(t)) u M u N dr, °T' MN =mj 8 D {x - V(r)) u' M u' N d T . (2.26) 

To calculate the leading order scalar bremsstrahlung it will be sufficient to know only the zeroth order term 
°h-MN of h,MN- So, in the sequel only °h,MN will appear and to simplify the notation, its left superscript will 
be omitted. 



2. The first order equations 

To derive the equations for the first order corrections to the particle world-lines one has to collect first order 
terms in the expansions of the embedding functions z M , z' M and the einbeins e, e' and choose suitable gauge 
condition to fix the r, r' reparametrization symmetries. From Eqs. (|2.9[) one finds for the first order corrections 
of the einbeins: 

l e = -co {x D ti MN u M u N + 2 , V = -e' (x D h M NU M u' N + 2 l z' M u' M ) . (2.27) 

The reparametrization freedom allows us to fix x e and V arbitrarily. We first substitute these expansions into 
Eq. (|2.8[) . collect all the first order terms and then choose the gauge fixing conditions x e = V = 0, that is 

x D h' MN u M u N + 2 x z M u M = (2.28) 

in the equation for to, and 

x D h M NU tM u' N + 2 %[U ,M = (2.29) 
in the equation for to'. The resulting equations for the first corrections to the particle trajectories read 

U MN %v = -x D Il MN (h' N r R h') u L u R , (2.30) 

















-xh,LR,N j 



R ,MN X-, n = _ X ^MN ^ nl r _ ) U ,L U ,R^ (2 31) 

where the projectors onto the space transverse to the world-lines are 

n MN = V MN _ U M U N^ jj/MN = ^MN _ JM^N ^ (2 32) 

whose presence demonstrates explicitly that only the transverse perturbations of the world-lines are physical. 



3. The second order equation for radiation 

The radiative component of the bulk scalar arises in the next order of iterations and is given by With 
appropriate combination of terms on the right hand side, Eqn. (|2.10[) is written as: 

U D ^{x, y) = j(x, y) = p(x, y) + a(x, y), (2.33) 

where the first term is localized on the world-line of the radiating particle to 

p(x, V) = ~f J d^(x -ur-b) S d (y)dT, (2.34) 

while the second is the non-local current 

a(x, y) = d M U' MN d N °$ - \h! d M °*\ , (2.35) 

with both h' MN and having support in the bulk. This current arises from non-linear terms in (|2.10[) due to 
the interaction of the bulk scalar with gravity. It has to be emphasized that it is non-zero outside the world-line 
not only on the brane but also in the bulk. Note that the decomposition into the local and non-local parts is 
ambiguous in the sense that part of the non-local term can be cast into a local form using the field equations. 
But the total source j never reduces to a local form as a whole. 
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D. The solution for °$, h MN and h M in M 4 x T d 

Our notation and conventions for Kaluza-Klein decomposition and Fourier transformation as applied to the 
ADD scenario with the transverse directions being circles with radii equal to R, are given in Appendix IA II 
It is important to stress at this point that in classical perturbation theory the interaction is described as an 
exchange of interaction modes (the classical analogs of virtual gravitons), but contrary to the Born amplitudes, 
where the simple pole diagrams diverge at high transverse momenta [3j , here the sum over these modes contains 
an intrinsic cut-off. Therefore the classical elastic scattering amplitude is finite and, furthermore, it reproduces 
the result of the eikonal method, if the eikonal is computed in the stationary phase approximation 25]. In 
other words, classical calculations in ADD correspond to non-perturbative ones in quantum theory (the eikonal 
method is equivalent to summation of the ladder diagrams). As will be explicitly demonstrated in the present 
paper, the same is true for bremsstrahlung. Specifically, it will be shown that the effective number of interaction 
modes is finite due to the cut-off and of order (R/b) d . Also, the summation over the emission modes is cut-off to 
a finite effective number of order (R'y/b)' 1 , which leads to a large extra enhancement factor 7 rf in ultrarelativistic 
collisions. 

Straightforward Fourier transform of (|2.22[) gives 

°$» = (2.36) 
P ~Pt 

Similarly, for the Hmn and H'mn, it is enough to Fourier transform the source terms (I2.26|) and plug into 
([2TTT) and (|2~T2"|) . One then obtains 

2Trx D mS(pu) l( b) f 1 \ 

KiNiP) = p 2_ p 2 T e [UMU N - jy—^miN J , (2.37) 

where p 2 — p^p^ , pu — p^u^ and p l T — ri 1 / R is the quantized momentum vector in the transverse directions. To 
get h /7 ^ N (p), one has to replace m — > m', u M — > u' M , b 0. Note that these fields do not describe radiation. 
They simply represent the scalar and gravitational potentials of the uniformly moving particles. Formally, 
this follows from the presence of the delta factors 8(pu) with pu — ^(p° — p z v ), from which it follows that 
p^p 1 ^ = p\iy 2 — 1) < 0, while the mass-shell condition for the emitted wave is p^p^ = p\ 5^ 0. 
Substitution of (|2.37p into (|2.30[) and integration of the resulting equation gives 

V*( T ) = y [ ^-JW e - f W (e-^ - l) (V M - t^« M - tJKp M ) , (2-38) 

w {2nfV^J (p 2 - p 2 T ){pu) \ Jy d + 2 2(puf )' y 1 

with 7^ = 7 2 — (d + 2) 1 and the Z?-dimensional vector p M = (p^,pj, = l l /R). It is easy to check that the 

gauge condition (|2.28p is satisfied. To ensure this exactly one has to keep the small second term in 7^, which, 

however, will be dropped in what follows in view of our interest in 7 > 1. We have chosen the initial value 

1 z M (0) = in order to preserve the meaning of & M as the impact parameter, namely 6 M = z M (0) — z' M (0). Note 

that the initial value of 1 i M (0) is non-zero and is computed from the gauge condition (|2.28j) . 

From (|2.38[) one can prove that the gravitational interaction does not expel the particles from the brane. 

Indeed, only the last j? M -term in the last parenthesis has non-zero components p l T = l l /R orthogonal to the 

brane. But the remaining expression is even under the reflection l l — > —l l and the sum vanishes giving z l = 
1 

The corresponding solution for the mass m' can be obtained by interchanging u M and u' M , replacing m by 
m! and omitting e _ ^ pb ' . 

E. <E> radiation - Basic formulae 

Finally, $— radiation is described by the wave equation (|2.33l) . which in terms of Kaluza-Klein modes is 

(□ + k 2 ) ^ n (x) = j n {x) = p n {x) + a n (x) , (2.39) 



1 This is on the average true also quantum mechanically. The <I>-quanta emission is symmetric on the average under reflection 
from the brane and the brane stays on the average at rest. However, to guarantee transverse momentum conservation in single 
4> emission in the bulk one should introduce brane position collective coordinates and deal also with the brane back reaction. 
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where k % T — n l /R, while p n (x) and a n (x) are 

p n {x) = -f J V(t) d^5 A {x -UT- b)dr (2.40) 

and 

= ^E^ (ftTW^^n-lM - ^'(X)^ . (2-41) 

respectively, with V 1 , °$ and ft.' given in (|2.38[) , (|2.36p and (|2.37l) . The rest of this paper is devoted to the 
solution of (|2.39|) and the analysis of the spectral and angular distribution of the emitted $— radiation. 

Once the solution of these equations is available, one can compute the energy and momentum radiated away 
using the standard formulae of radiation theory. To compute the momentum loss due to scalar bremsstrahlung 
emitted during the collision one considers the world tube with topology R 1 ' 3 x T d , with boundary dVt = 
EooUS-ooU-B consisting of two space-like hypersurfaces S±oo in R 1 ' 3 &tt = ±oo and the time-like hypersurface 
B at r — > oo and integrate the difference of the fluxes through E±oo to obtain 

P^ = Jd d y( J - J \T v »d 3 Y, v . (2.42) 
V £_«/ 

Here one makes use of the brane components of the energy-momentum tensor of the bulk scalar (it is easy to 
show that there is no radiation flux into the compact dimensions |23[). Start with 

t mn = qM^qNq _ l f? jifjv (e$ ja i (2 43) 

where only ^ has to be taken into account, since °$ is not related to radiation. The integral over B is zero 
due to fall-off conditions, so the difference of the surface integrals (|2.42p can be transformed by Gauss' theorem 
to the volume integral 2 

P» = J d d y J d N T N »d 4 x = Jd d y J D D $ d 4 x. (2.44) 

v n v n 

Using f|2 .33[) and the retarded Green's function of the D— dimensional D'Alembert operator to solve the wave 
equation 

one obtains 

d 3 k 



fe0 = ^k 2 +fc| 

where j n (k) is the Fourier-transform of j n (x) (for precise definition see Appendix I A II) and will be referred to 
as the radiation amplitude. With the parametrization k = |k| (sin cost/?, sin#sini/2, cos#) one obtains for the 
spectral-angular distribution of the emitted energy E — P°: 



dE 



^— Vk 2 |j n (fc)| 2 , dn 2 =sm9ddd<p. (2.47) 



d\k\dn 2 16tt 3 V 

Therefore, the leading order radiation loss is determined by the Fourier-transform j n (k) of the source in the 
four-dimensional wave equation (|2 .33[) for 1 <£>" on the mass shell of emitted waves 

= k 2 T . (2.48) 

If the impact parameter is small compared to the compactification radius (b <C R), the summation over KK 
masses can be replaced by integration according to (|A.9[) . with integration measure d d kr = k^T 1 dkx dQd-i- 



2 It was taken into account that T d has no boundary 
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Figure 1: The angles in lab frame used in the text. 



Since the radiation amplitude j n depends only on |fex|, integration over the angles is trivial and gives the volume 
Qd-i of the unit d — 1— dimensional sphere. Therefore, 



dE n d 



d- 1 / l,2,rf-l 



d\k\dn 2 2(2tt) 



D-l / T 



k 2 k!f L \j n (k)\ 2 dk T . (2.49) 



with n in j n expressed in terms of kr, and k° = ui = y'k 2 + k'i,. 

When the summation over KK emission modes is replaced by integration, one can compute the emitted 
energy using directly higher dimensional Minkowski coordinates. Take the coordinate system with angles on 
the (D — 2)— dimensional sphere Q.d_2, with i? the angle between the (D — 1)— dimensional vectors K = (k, k l T ) 
and u, <j) the polar angle in the plane perpendicular to u varying from (direction of b) to 2ir (see Figure []} 
and integrate over the angles to obtain |23| 

dE UJ^ 2 . , , , ,n , 

\j(k)\ 2 . (2.50) 



dwdQ d+2 2(2tt 



irf+3 1 



where u = (ku') = k° — \J\s? + kj, is the higher-dimensional frequency. 



3. THE RADIATION AMPLITUDE 



As we have seen, the radiation amplitude j n (k) is the sum of a local p n (k) and a non-local part a n {k). 
These two parts have different dependence on the frequency and the angle 9 of the emitted wave with respect 
to the direction of collision. Both have an intrinsic cut-off at some (angle-dependent) frequency, which in the 
ultrarelativistic case is high compared to the inverse impact parameter 1/6. Typically, the two amplitudes cancel 
each other in some range of angles and frequencies. To obtain the correct expression for their sum one has to 
carefully take into account non-leading contributions. 



A. The local amplitude 



Fourier-transformation of (|2.40[) gives 



P n (k)=if 



— i f p i( kb ) 



e l(fctl)T (fc 1 z)dr 



(3.1) 



where the scalar products are denoted as (ab) = a^b 11 and for simplicity the parentheses will also be omitted 
when it is not ambiguous. Substitution of ^ from (|2.38|) and integration over r gives 3 : 



P n (k) 



x- ri m'fv nh!A 



A-K 2 V{kuf 



E 



ku ( jku 



ku 
d + 2 



-k 



(3.2) 



3 Note that the unity inside the first parenthesis of I I2.38H corresponds to a constant 1 z and does not contribute to radiation. 
Formally, its contribution to p n vanishes by symmetric integration 12311 - 
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where the integrals /; and If are defined by 



S(pu') 8{ku - pu) e- i( - pb 1 4 _ 
a p, L 



p 2 -p\ 



6(pu') 6(ku - pu) e~ 1 ^ 



P 2 -Pt 



P 11 d A p, 



(3.3) 



The sum over / represents the sum over the interaction modes labeled by the set of integers l l , while the 
dependence of the amplitude p n {k) on the vector index n, which labels the emission modes n l , is hidden inside 
the fc°— component of the wave vector: fc° = ^k 2 + k\. The integrals are given in [23| and lead to Macdonald 
functions: 



2tt 



with 



Ko(zt), If 



(ku)b 



2tt 
^vb 2 

_ (ku')b 



bzKo(zi) 



+ 1X^1)^ , 



, Zl ^{z 2 +p 2 T b 2 f' 2 



(3.4) 



(3.5) 



and the hatted Macdonald functions defined by K v {x) = x v K v (x) and having for i/^O a finite non-zero limit 
as x — > 0. Thus, in terms of Macdonald functions the local amplitude is: 



P n (k) = 



x 2 D m'f 



AttvV 



i(kb) 



E 



2- 



? 2 

v '7 



1 



1Z 



j\d + 2 2v 2 j 2 



K (zi) - i 



-ll (kb) 



2 2 2 

7 Z v jz 



K x { Zl ) 



(3.6) 



1. The 7 — > oo limit. Mode and frequency cut-offs 
In the ultrarelativistic limit 7 — > 00 the leading terms of p n are 



*D™'f „i(kb)X^ \ z 'ts r.A A kb ) t><„a d+1 f z ' T d \ t.~ 1 ■ 1 ...o,.-^-- (3.7) 



This is a systematic ultra-relativistic expansion in powers of I/7, modulo the coefficients of the various Mac- 
donalds as well as the overall factor in front, which depend on the velocity v. However, as will become evident 
below, this form is adequate for the following discussion and the computation of the emitted energy to leading 
order. 

Two important general remarks are in order here: 

(a) The effective number iVj nt of interaction modes. The exponential fall-off of the Macdonald functions at 
large values of the argument zi leads to an effective cut-off ./Vi nt of the number of interaction modes I in the sum. 
One can estimate the radius p™* of the sphere in the space of l % , beyond which the modes can be neglected, by 
setting 

(p^) 2 6 2 ~l, (3.8) 

from which 

R\ d V 



N i*~\j) ~I3- (3-9) 



For iVi n t 3> 1, which is the case of interest here, one may use (jX9|) to obtain (Z > 0) 

1 £ K X (y/Z»+ft?) - ^j^K x+d/2 (Z) (3.10) 



and end up with 



P n (k) - -^P- {7 £./*(*) i^K d/2+1 {z) + (d ^±±^) K d/2 (z) + "0(7- 2 )'j , (3.11) 



where 



x 2 D m'f 



2{2it) d / 2 + 1 b d 



(3.12) 
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(b) Angular and frequency characteristics. The local radiation amplitude above in the b <C R limit is expressed 
solely in terms of Macdonald functions with argument z. Later, it will be shown that the non-local amplitude 
contains also Macdonald functions but with argument z' . The exponential fall-off of these functions implies the 
effective cut-offs z ~ 1 and z' ~ 1 in the corresponding radiation amplitudes. These, in turn, translate into 
angular and frequency characteristics of the corresponding radiation. 

Specifically, with 9, a and $ as shown in Fig. [IJ define 

ifj = 1 — v cos 9 cos a = 1 — v cos$ , (3.13) 

which satisfies 

J =li>, (3.14) 

and in the ultrarelativistic limit varies in the interval I/27 2 ~ 1 — v^ip^l + v~2. 

Consider the neighborhood of z *~ 1 which gives the dominant contribution of the local radiation amplitude. 
One has to distinguish various domains of the emission angles. For small emission angles 9, a one has 

1>~ lh' 2 + 2 + a 2 ), (3.15) 

so that (i) inside the small cone 9 2 + a 2 < I/7 2 one obtains if) ~ I/7 2 , so that the characteristic frequencies 
oj ~ 7 2 jb. This angle- frequency regime will be called in the sequel the z-region. (ii) For if) ~ 1, i.e. for 
a, 9 ~ C(l), one obtains z' ~ z/7, which implies a low frequency regime uj ~ 1/6, whose contribution to the 
emitted energy is negligible in view of the relative smallness of the phase-space factor in (|2.49[) . 

To summarize, the above analysis of p n combined with the phase space factors in (|2.49p . leads to the conclusion 
that the leading contribution of the radiation due to the local amplitude is beamed, i.e. directed inside the small 
cone 9 2 + a 2 < I/7 2 and has high frequencies u> ~ J 2 /b. Radiation with these characteristics will occasionally 
be called z-type. 



B. The non-local amplitude 

The non-local amplitude obtained from (|2.4ip by Fourier transform is 

J2- ™' ffL„J\2 



with 



^fc) = ^e«J»W, J«(fc) = ^J2j nl (k), (3.16) 



{k) -J P (p 2 ~P 2 T )[(k- P ) 2 -(k T -p T ) 2 }- ^ 



k° = \/M 2 + k\ and k is a 3— dimensional vector lying on the 3— brane, where k l T = n l /R and p l T — l l / R with 
integers {n l },{l 1 } are d— dimensional discrete vectors corresponding to the emission and interaction modes, 
respectively. 

Using Feynman parametrization J nl takes the form: 



1 

jnl 



J [p 2 - (k T x - p T ) 2 } 2 

u 

Integrating over p° and splitting p into the longitudinal pn and transversal pj^ parts, integrate over pii. Then, 
introducing in pj_ the spherical coordinates c? 2 pj_ = |p±M^irf|pj_| and integrating first over the angles and then 
over |pj_|, one obtains 



1 

jnl = ?L f dxe^ x K^(Cni), 

IV J 



(3.18) 



with 

C 2 i(a;) = z' 2 x 2 + 2 7 zz'x(l - x) + z 2 {\ - x) 2 + b 2 {k T x - p T ) 2 . (3.19) 
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Again, the summation over I is performed trivially for b -C R using (|3.10[) 4 . The result is 

1 

J n (k) = A d J dxe-'W'Ky^ i(Cn); A d = (27r)d/2 ■ (3.20) 
o 

with 

£(x) = z' 2 x 2 + 2 1 zz'x(l-x)+z 2 {l-x) 2 ; ( n (0) = z , ( n (l) = z' . (3.21) 
Writing £ 2 successively in the form 

C 2 (x) = -fx 2 + 2/3x + z 2 = a 2 -r 2 , (3.22) 

with 

£ 2 = 2 7 zz' - z 2 - z' 2 = lu 2 b 2 sin 2 (9 cos 2 a + b 2 k 2 T = z /2 7 V sin 2 ■& , /? = 722' - z 2 , (3.23) 

and 



2 + z^ , r = i ( x - £ ) , (3.24) 



and using formula f. 2. 16. 12-4 of [36[ 

2 l/2 



Kv-i/2\yJa 2 -r 2 \ = - T j 1 a 2v I cosh(ry)K- v ( a ^ y 2 + lj dy , o -1 and a > (3.25) 

o 

for jtx =— 1/2, v = (d — l)/2 one may rewrite (|3.20[) in the form 

1/2 2 °° ^ 

j»(fc) = A d ^- T ^ /" dyk- u {a^ftF+l) J dxe-^ kb >cosh(ry). (3.26) 
o o 

Perform, next, the integration over x and introduce the additional angle </>, so that the generic (D — 1)- 
dimensional unit vector K/|K| (the normalized higher-dimensional emission vector K) is decomposed as: 

-^7 = 7^- cos $ + T-^T sin d cos d> + m sin # sin 0, (3.27) 
|K| |u| |b| 

where m is a D — 1 dimensional unit vector orthogonal to the collision plane (spanned by u and b). Then 
(k ■ b) — — 7z't> sin t9 cos = — £cos</> and a = ubip/ sin-d. Substituting this into (|3.26l) we have 



ifi/i f ^ — t ] \ 1 y z + cos z 

s 3=0,1 I y y 

(3.28) 



where Sj = j — (3/S, 2 , j — 0, 1. 5 The convergence of these integrals is controlled by the competition of the 
exponential decay of the Macdonald function and the exponential growth of the hyperbolic functions. In all 
cases the first is faster, but when the difference of the two arguments is small, the main contribution to the 
integral over y comes from large values of y. 



4 Using (IA.9I ) the summation is converted to integration over d d l. One then shifts the integration variable p' T = pj< — xkj< and 
applies 113.101 1. In the present case Z is not constant but depends on x. However, as it can also be checked numerically, 113.1011 
and the subsequent treatment is a good approximation for any x ^ 1, because Z(x) > 1 in all relevant frequency regimes. 

i 

5 Direct integration of l3~20l) for (0 = 0, k T = 0) gives / K i/2 _ 1 ( ^2fix + z 2 ) dx = P' 1 [K d/2 {z) - K d/2 (*')]• The constants of 



integration of the terms j = 0,1 are chosen so that for 9 = they satisfy J"|(8 = o k T= o) = A,j/3 — 1 -K*d/2( 2 )i J£/ \(e=0,h T =0) = 
-A d ^k d/2 (z>). 
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Since y 2 + 1 1 and ^ sin 2 ^ 1, one can equivalently write 



j;(A0= (-iJi+ie-yWA,,-^ sin 2 V y rfy 2 V + 1)t//2+1 / fc [ysiah(^-y) - i cos cosh 



The y-integration for each value of k is performed by successive applications of the identity 34 [ 



K v (z) = K v+2 {z) - 2{V + 1) K v+1 {z) (3.29) 



in combination with 

a v Idy V 



(S +1) I COShtfJj!,) I ° I I 



obtained from (|3.25|) . with argument Zj = y a 2 — ^ 2 <5? , i.e. zo — z, Z\ 
For example, using p. 291) the k = term leads to the integrals 



a v \ dy 



y+2 



( y2 + 1 ^/2+l a ( y 2 + 1)^+3)72 



y sinh^y) 
cosh^y) 



(3.31) 



The first term in the square brackets is given by (|3.30[) . The second is computed using again (|3. 291) . which leads 
to two new integrals, the first of which is 

J K v+3 (a^TT) j y ^ H(6 . y) 1 i J ^^+3/2^) 1 

a (y2 + 1)(,+ 3 )/2 s i a 4 i * r n r 

suppressed for a 3> 1 compared to ()3.30j) . Similarly, the second is further suppressed by two more powers of a. 
Terms with increasing k are evaluated in the same way and lead to further suppression by inverse powers of a 2 . 
The end result for Jo(k) keeping terms up to 1/a 4 is then 

W = Jfr (pKdM*) ~ i(kb)K d/2+1 (z) - ^±plR d/2+1 (z)+ ^±K d/2+2 {z)^ + R z . (3.33) 

Notice that J${k) is a series of Macdonalds with argument z. Consequently, it is important mainly in the 
z-region, where a = wbip/ smfl ~7> 1, a self-consistency check of our approximations. The coefficients of all 
terms in (|3.33[) have expansions in powers of 7~ 4 . In the z-region the first term starts with 0(7~ 3 ), the second 
with 0(7 -4 ), the third and fourth terms with order 0(j~ 5 ), while the remainder R z = (D(-f~ 6 ). 

Following the same procedure J™(fc) is written as a series of Macdonald functions with argument z', namely 

Jf(fe) ~ A d e-^ ft K d/2 (z>) - K d/2+ i(z')) + R z > , (3.34) 

whose main contribution comes from the region with z' ~ 1, i.e. u> ~ 7/6, # ~ 1, in which indeed a~7> 1. 

The condition a 2 ^> 1 is not satisfied in the region with 8 ~ However, in that region both the exact 

expression and the approximate one have negligible contribution to the amplitude. Figure [2] displays graphically 
the maximal difference in the real part of J" (k) between the two expressions. 

Incidentally, notice that there is no strong anisotropy in cj> in the z'— region: the real part of main terms of J" 
(|3.34[) is independent on 0, while its imaginar y p art depends only by the overall factor cos <j). The same picture 
was obtained without any approximations in [23j, where scalar mediated collisions were studied. 

Going back to (|3.16|) one sees that <j n (k) is the sum of two sets of Macdonald functions, one with argument z 
and the other with argument z' . In analogy with p n , the first sum contributes mainly in the z-region. Similarly, 
the leading contribution of the second set of Macdonalds comes from the region with z' = wb/^fv ~ 0(1). In 
the ultrarelativistic limit this translates into angular and frequency characteristics of the emitted radiation. For 
generic values of the angles, this means u> ~ 7/6 and defines what we will call the z'— region and, correspondingly, 
z'—type radiation. It is unbeamed radiation (t? ~ 1) with characteristic frequency uj ~ 7/6. 

The non-local pieces (fc) and cr™(fc). It is convenient to separate the two kinds of contributions to the 
non-local amplitude by writing 

<T n (k) = a%(k) + o?(k), (3.35) 
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Figure 2: The real part of the amplitude J™ (solid line) in D = 4 dimensions, calculated numerically from the exact 
formula (|3.28f) . compared with the approximation given by (|3.34p (dashed) for 7 = 1000 and for cos(j) — 1, which gives 
the largest deviation between the two. For 7# > 4 the deviation is negligible. For 7$ < 4 the exact expression will be 
used numerically. The imaginary part has similar behavior, and furthermore is suppressed compared to the real part by 
inverse powers of 7. 



with the first (second) given by (13. 16)) with Jfi (J") on the right hand side. Thus, 

o$(k) = Ae'W Tj^ (pk d/2 {z) - i(kb)K d/2+1 (z) - il±iM/? d/2+1 (z)+ ^L±K d/2+2 (z)^ , (3.36) 

and 

<(fc) ~ [{e - !3) K d/2 (z') + i(kb) K d/2+1 {z')) , (3.37) 

respectively. 

Correspondingly, the total radiation amplitude j n (k) is written as a sum of two terms, one function of z, and 
the other function of z' 

r(k)=j:(k)+&(k), j2(k)=p n (k)+aZ(k), rt(k) = a[ l (k) (3.38) 



C. The part j"(k) of the radiation amplitude and destructive interference 

1. j™ in the frequency range w 2> 7/b 



Consider first the regime with 1? ~ 1. Here z ~ 7 and from (|3.36[) and (|3.1ip one obtains that j™ ~ exp(— 7) 
due to the Macdonald functions. 

Now take the most interesting case of ■& ~ I/7, in which z ~ 1. Add p™(fc) and (Jg (fc) given in p. lip and 
(|3.36p . respectively, and use the ultra- relativistic expansions to obtain in leading order: 



j n z {k) 



A(d+l)e 



i(kb) 



d + 2 



k d/2 (z) 



cos 2 a 



(sin 2 6> + tan 2 a) K, 



sin 2 sin 2 <p 



d/2+1 



(2) 



1 



A", 



2/ 2 + 2 



(3.39) 



All terms inside the square brackets are of 0(7 -2 ). Given that in the z-region l/jip = z '/ z ~ C(t)> the 
leading contribution to j™ above is of C(7 _1 ). Higher order terms have been ignored. The terms of order 0(7) 
and O(l), both present in the ultra-relativistic expansions of <7q (fc) and p n (k), have opposite signs and cancel 
in the sum. This is a general phenomenon of destructive interference related to the gravitational interaction. 
Thus, the two leading powers in the ultra-relativistic expansion of the direct $— emission amplitude from the 
accelerated charged particle, cancel the ones coming from the indirect emission due to the $ — $ — h interaction. 
As a consequence, the z-type (beamed and high frequency) part of the radiation is highly suppressed in the 
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ultra-relativistic limit, compared to the naive expectation. One can check that destructive interference is valid 
also in the case of $— radiation in arbitrary D— dimensional Minkowski space-time, which can be obtained as a 
limit of the present discussion. It was first observed for gravitational radiation in D — 4 [3(| (using a different 
approach) and it was recently generalized to arbitrary dimensions in [24[ . For the system at hand, an alternative 
proof is presented in Appendix El using a different approach also suitable to the frequency range ui ^> j/b. 

The following comments are in order here: (a) As a check of the above series of approximations, one may 
consider the special case of 6 = 0, for which (kb) = — £ = 0. In this case the exact value of J™ obtained from 
(|3.20[) coincides with the one obtained from the approximate expressions (|3.33l) and (I3.34[) . (b) Furthermore, 
equation (|3.39[) can be shown to coincide with the corresponding quantity in the case of non-compactificd 
D = A+d— dimensional Minkowski space. This generalizes to scattering and radiation processes in the relativistic 
case, the non-relativistic argument about the behavior of Newton's potential, i.e. that at distances b <C R a 
pomt particle generates the D— dimensional potential, while at b 3> R its potential behaves as four-dimensional 

One may equivalently parametrize using the angles $ and <f> and write it in the form: 

1 ,~ , , , sin 2 tf 



\ i(kb) 



d + 2 



sin u f *■ \ 

(2tf - 7 - 2 )^ /2 (z) - — — [{d + l)K d/2+1 (z) - sin 2 (f>K d/2+2 (z) J 



(3.40) 



Note that in the computation of the emitted energy below both angles will be taken continuous; a sensible 
approximation for N- lriii 3> 1 assumed here. 



2. j" in the frequency range lo < ■y/b 

For uj < 7/6 and d ~ l/ 7 using (|3TTTT) for p n and (TXToT) and (j3~2U|) for a n , one concludes that \p n \ > \cr r ' 
and, therefore, 



J n (k) 



P n {k) 



1 y . sin § COS0 y 

^ Kd ^ z)+l -^b~ Kd ^ z) 



(3.41) 



For § ru 1, on the other hand, p n , Cg and <r™ are all of the same order, but suppressed compared to the previous 
case. In addition, the contribution of this regime to the emitted energy will be shown to be further suppressed 
by the integration measure. 

More interesting is the case with cj ~ 7/6. If # ~ 1, then z ~ 7 and using (I3.11|) and (|3.36|) one concludes 
that j™ is exponentially suppressed because of the Macdonald functions. However, for i3 ~ 1/7, one may use 
([3~TTj) and (pHSl) to obtain that p n ~ 7 and ct# ~ 7, respectively 6 . 



D. The part of the amplitude 



Equation (|3~37)) can equivalently be written in the form: 



1 \ ■ ... . sin 9 cos a cos ip ~ . 
-ay - 1 K d/2 (z') + 1 ^K d/2+1 (z') 



Furthermore, using the angles 1) and it becomes: 



Jz 



A 



1 



\7 w 



l)K d/2 (z') + i- 



n $ cosq 
'yz'ip 



-K, 



d/2- 



-I CO 



(3.42) 



(3.43) 



First, for cu ^> 7/6 one has z 1 ^> 1 and, consequently, j™, in (|3. 431) is exponentially suppressed. Next, for 
(cj ~ 7/6, ■& ~ 1) j" in (|3.43p is dominated by its real part which is of order 0(1/7). For (u <C 7/6, i9 ~ 1) one 
obtains j™, ~ <7q ~ 1/7. As will be shown below, however, this region contributes negligibly little to the emitted 
energy. Similarly, for (w <§; 7/6, ^ I/7) on the basis of (|3.11|) and (|3.20p one concludes that the amplitude 
j™ ~ ctq <C p" ~ 7 2 . Finally, based on numerical study and previous results in D = 4 [3(j one obtains that 
(|3.43|) is valid also in the regime (a; ~ -f/b, ~ 1/7) and gives j™ ~ 7. 



Using the formulae of Appendix IA 31 one gets in this kinematical regime a = z'^yip/ sin $ ~ 1 and also £ ~ 1 as well as /3 ~ 1. The 
integrand in 113. 28H is independent of 7. All 7 dependence comes from the overall coefficients in I I3.28H and 113. 1611 . 
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Figure 3: Frequency (a) and angular (b) distribution for d = and 7 = 10 5 



E. Summary 

The behavior of the local and non-local currents in all characteristic frequency and angular regimes is sum- 
marized in the following Table I. 



u> <C 7/6 



uj ~ 7/6 



u > 7/6 



no destructive interference 
j n ~ p n > a% ~ trf 



no destructive interference 
# ~ p" ~ a " ~ 7 

j™ = (333) ~ 7 



destructive interference 

# - eaa ~ p"/ 7 2 ~ 1/7 

# - (Eli ~ exp(- 7 ) 



no destructive interference 

,'D _j n n _j „-n 



destructive interference 

# - (EU+rai) ~ exp(- 7 ) 



7 



destructive interference 

i" - mn^nm ~ e XP (- 7 ) 

# = (EH ~ exp(- 7 ) 



THE EMITTED ENERGY - SPECTRAL AND ANGULAR DISTRIBUTION 



The spectral and angular distribution of the emitted energy is obtained from (I2.49[) or (|2.50[) . The integrand is 
the sum of three pieces proportional to | j™ (k) | 2 , | j™ (k) | 2 and j ™ j™, +jzj™ (the bar denotes complex conjugation) , 
so the total radiated energy splits into three parts 



dE = dE z 



dE z 



dE z 



(4.1) 



which will be called z— , z' — and zz'— radiation, respectively. All terms contain a factor of A 2 from the square 
of the amplitude, an explicit l/2(2-7r)( d+3 ' from the Fourier transform in (|2.50p . and a factor & _ ( d + 3 ) from the 
corresponding power of u in the integrand. This leads to a general expression for the total emitted energy of 
the form 



E 



1 



4 /2 rO 



8 ( 27r )2d+5 6 3d+3 



■T 



(4.2) 



with an overall coefficient expected to be of order one and the power of gamma depending on the particular 
type of radiation under discussion and which is easily determined as follows: As argued above and in [3(| 
and also shown for example in Figure |3] (obtained numerically), the $— radiation is emitted predominantly in 
well-defined relatively narrow frequency and angular windows, and with amplitudes shown in Table I. Thus, it 
is straightforward to estimate the powers of 7 in the various components of its energy, using (|2.50[) 



E 



duj / <M|j'| 2 w 



d+2 ■ d+1 , 



(4.3) 
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with \j\ 2 being |j™| 2 or | 2 or + and with the range of integration not contributing extra factors 

of 7. For example, the contribution of j", which is dominant in the regime (u) ~ 7 2 /&, 1? ~ I/7) has I/7 2 from 
|j"| 2 , (l/7) d+2 from the angular integration, and (7 2 ) d + 3 from the integration over uj, with the final estimate 
being r y d+2 . 

The result of this computation is the content of Table II below 7 . 



\ UJ 
V\ 


uj <C 7/6 


UJ ~ 7/6 


cj ~ 7 2 /6 


UJ > 7 2 /6 




negligible 
(phase space) 


Ed ~ "f 3 , from j" and j™, 


£ d ~7 d+2 , fromj™ 


negligible radiation 


1 


negligible 
(phase space) 


£<i ~ 7 d+ \ from 


negligible radiation 


negligible radiation 



We proceed next to the detailed study of the various components of radiation with the frequency and angular 
characteristics of the three most important cells of Table II. 

A. The z— type component of radiation with uj ~ j 2 /b 

According to Table II, the z-type radiation (due to |j"| 2 ) is always beamed inside •& ~ I/7. Furthermore, 
for d ^ 2 it is dominant with characteristic frequency uj ~ 7 2 /6. The cases d = and d — 1 will be treated 
separately in another subsection. 

It is convenient to write the current j z (|3. 401) in the form 



s=0 

A d+ 1 



J 7^ + 2 
V* = -A (d + 1 



1 

7 2, 



sin & 

~fipz 2 



K d /2+i{z) 



Jz 



sin 2 tfsin 2 . 

A — — K d /2+2\z) 



Squaring and substituting into (|2.50[) one obtains 

dE z 
dujdQd+ 



,d+2 



2(2tt) 



d+3 



E a 

a,b=Q 



Jz jz 



To integrate over frequencies it is convenient to change variable from uj to z and define the quantities 

C$ = f kd /2 ^z)kd, 2+b {z)z d+2 ^ +s ^dz. 



(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 

(4.9) 



The integration over 1) is performed using (|A.11[) . Finally, the integration of s'm 2 <p and sin <fi over the remaining 
angles of S d+1 is f^+i/2 and 3f2<2+i/8, respectively. The angular and frequency profiles of this component of 
radiation in dimensions d 2, for which it is dominant, were obtained analytically and numerically, respectively, 
and have the general form shown for d = 3 in Figure [4] 

The end result for the total emitted energy in this component of radiation is 



E z 



2(27r6) d + 3 



^ 2^ °a6 U ab 



(4.10) 



.,6=0 



7 Note that in D = 4 it is found that all three types of radiation are of equal 0(-f 3 ). This seems to disagree with ;30], where it is 
stated that the leading C(7 3 ) is due to z'— type alone. 
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dE z 
y 2 d(0 



(ob/y 




(a) 



(b) 



Figure 4: Frequency (a) and angular (b) distribution of z-radiation for d = 3 and 7 = 10 3 



where 



D 



(d) 
00 



2 d+1 (d+l) 2 r 2 (^) 



D 



T{d + A) 

_ 2^+ 4 (rf + i) 2 r(^)r(^) 



(d) 
01 



2 d+3 (d + i) 2 r 2 (^) 

(d + 2)r(d + 4) 



-D, 



(d) 

02 



11 



r(d + 5) 



(d) 

12 



2(d+l) 



(d) 

11 > 
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f £,W 

2(d+l) 01 ' 

(d) 



8(d + l) 



2^11 



(4.11) 



When integrating over z in (|4.9[) one should remember that the expansion (I3.40[) is accurate in the high 
frequency domain, around and beyond z ~ 1. However, for rf ^ 2 it can be checked both analytically and 
numerically that the integral from to 1 of the difference of the exact energy density based on (|3.11l) and (|3.28p 
and the approximate one based on (|4.4[) is negligible. Thus, one can conveniently expand the integration region 

in (|4.9p from to 00 and evaluate using (|A.13I) . 

Collecting all contributions one obtains for the energy of high frequency z— type radiation 



_ r < >4> m ' 2 / 2 „,d+2 
E - Cd b 3d+3 7 

with C 2 = 1.42 x 1(T 6 , C 3 = 6.02 x 1(T 7 , C 4 = 3.45 x 1CT 7 , C 5 = 2.67 x 1CT 7 and C 6 = 2.76 x 1(T 7 . 



(4.12) 



B. The z'— type radiation with # ~ 1 



According to Table II, wide angle radiation (1? ~ 1) is mainly z'-type (due to |j™/| 2 ) in all dimensions and 
has characteristic frequency oj ~ 7/6. Also, for d 3 radiation with u ~ 7/6 is predominantly emitted in wide 
angles. 

Squaring (I3.43[) , substituting into (12.50[) and integrating over w from to 00 and all angles except 1?, one gets 
the angular distribution 



dE*' _ ^m' 2 jy+ 1 r(^3)r 2 (^)r(^) S in d+1 ^ 
M - b 3d+a 2 7 7r 3d / 2 + 4 r(^)r(2d + 3) V 2 



(4.13) 



Formula (|4. 1 3[) gives the dominant wide angle radiation in all dimensions d J? 3. Figure [5] shows the angular 
and frequency profile of this component of radiation for d = 3. To compute the total energy of this type we 
integrate over 1? making use of (IA.12|) . 



For d ^ 3 the emitted energy is given by 



E* 



, K D m! f d+1 



c', 



2 d-8 T (3d±3) r 2 (^) r(^) r (^) 
7r 3d / 2 + 4 r(2d + 3)r(d) 



For d = 2 one obtains 



E z 



105>f 4 m' 2 / 2 - 
— 7 J In 7 

2!6(27T)769 ' ' 



(4.14) 



(4.15) 



The cases d = 0, 1 have to be considered separately since for them z-, z'- and zz'-types of radiation are 
comparable and splitting the amplitude into j z and j z > is not particularly useful. 
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Figure 6: (a) The ^-distribution in 4D for 7 = 10 4 and (b) in 6D for 7 = 10 3 . 



C. The cases d=0, 1 



According to Table II the emitted energy in 4D is concentrated in the region to ~ 7/6, 6 ~ 7 _1 . In this case 
the exponent e 1 x ub sm * cos v in the stress amplitude cr(fc) does not oscillate fast and the emitted energy may be 
easily computed numerically. The frequency and ^-distributions in this case are shown for 7 = 10 5 in Figure 
while the distribution over <j) (which coincides with tp in 4D) is presented by Figure [6] 

The total emitted energy is 

E = C * 4 ™ 1 7 3 , C «8.3xl0- 5 . (4.16) 
The frequency distribution is non-zero at to = (see Figured]), in agreement with the analytically derived value 



1 xfm' 2 f 2 2 

7 7 



w=0 3 x 2 6 x 7T 4 6 2 



due mainly to the imaginary part of the p— amplitude. 

The frequency distribution of the emitted energy E\ in 5D for 7 = 10 4 is shown in Figure[7] It is characterized 
by a long tail beyond the value uj ~ 7/6, which as can also be argued analytically 8 leads to a behavior 
dEi/duj ~ 1/lo (Figure [7(b) \ all the way to uj ~ 7 2 /&, beyond which it falls-off exponentially. The integral 



Notice from Table I that the total amplitude satisfies j(ui ~ 7/6) ~ 7 2 j(kJ ~ 7 2 /6). This gives for ci = 1 the estimate 

\j\ 2 LJ d + 2 ~ 
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JE_ 
dm 



dK_ 
dm 



cob/ 7 



(a) 



cob 



(b) Intermediate region uib = 2O7 — O.I7 2 : exact distribution 
(solid black) versus the approximate dE/dui ~ l/w (dashed 
red) 



Figure 7: Frequency distribution for d = 1 and 7 = 10 4 : (a) for uj < 7/6, (b) for 7/6 < lo < 7 2 /6 



of dEi/du> over the range (7/6, 7 2 /6) gives an extra logarithm in the total emitted energy, which is computed 
numerically to be 



4 /2 £ 2 

E x = d * 5 ™ 1 7 3 ln 7 , Ci w 1.64 x 10" 5 . 
0° 



(4.17) 



D. The estimate of the zz' — interference part of radiation 



The purpose of this subsection is to estimate the contribution of the interference part (j™ j™, + c.c). It will 
be shown that it is subleading for d ^ 2 and of the same order as z— and z'— contributions for d = 0, 1. 

The interference term E zz ~ /0™ij + c.c.) w rf+2 <ir2d+2 contains the product of Macdonald functions 
i^(z)-ftT(z'). Thus, its value depends on the overlap of these functions in the domain (z < 1, z' < 1), or 
equivalently (to < 7/6, < 1/^/7). The presence of the factor u d+2 implies that most of the contribution to 
the integral comes from the large uj regime with w ~ 7/6, in which z 1 ~ 1. 

For z <C I/7 the integral is suppressed by the volume factor. Thus, the interesting regime of z is 7 _1 < z < 1. 
In this regime one may estimate the contribution to the interference integral using 



Xe i{kb) gin 2^ 

2 j?pz 



2 K d/2+l 



(z). 



3z 



X 

1^' 



K d/2 (z'), 



(4.18) 



since, it can be checked from (I3.11|) . p. 361) and p.43p . that they are either dominant or of the same order as 
the remaining terms. 

One then obtains for the interference part of the energy loss 



dE z 



A 2 cos(w6 sin d cos < 



Integration over all angles except i9 gives 



d+3^2^2^2 



-K d/2+1 (z)K d/2 (z') sin 2 tfw 



d+2 



dE z 



A 2 



dti duo 2 d + 2 7r( d + 3 )/ 2 r((d + 1)/2) 7 V V 



■ J (u b smfl) K d/2+1 (z) K d/2 (z') uj 



d+2 sin d+3 { 



(4.19) 



(4.20) 



The value of the integral over 1) and w is controlled by Jo- For i9 ~ l/-y/7; z an d z ' are both of 0(1), while the 
argument of Jo is of 0{ y / 7 ~f) 3> 1. Using then the asymptotic expansion of Jo and approximating the hatted 

Macdonalds by their values at z ~ z' ~ 1, one can estimate as in previous cases the power of 7 in E zz to be 
7 d / 2+3 / 4 . This is negligible, compared to the other contributions in all dimensions. 

One is left with the contribution from d ~ I/7, where ui b sin?? ~ 1 and z ~ 1/7. Substituting Jo ~ Jo(0) ~ 1 
and K(z) ~ K(z — 0) one estimates the integral to be of C(7 3 ) in all dimensions. This is of the same dominant 
order as E z and E z for d — and d — 1 and was included in the total energy evaluated in the previous 
subsection. 
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5. SUMMARY OF RESULTS 

Scalar bremsstrahlung radiation during the transplanckian collision of two gravitating massive point particles 
in arbitrary dimensions was studied classically in the laboratory frame. The main goal was to compute the 
powers of the Lorentz factor 7 and how they depend on the number of extra dimensions d. 

We computed both analytically and numerically radiation into truly massless and massive (for the branc 
observer) modes. An essential difference with the previously considered case of scalar bremsstrahlung in flat 
space from particles interacting via a scalar field, is that in the latter the scalar field(s) is linear, while here the 
bulk scalar interacts non-linearly with gravity. Within the perturbation theory with respect to both scalar and 
gravitational coupling constants it was found that the radiation amplitude consists of a local and a non-local 
part. Furthermore, it was shown that in a certain range of angles and frequencies the leading terms of these two 
mutually cancel, while the remaining terms can be presented as the sum of two contributions (j z , j z >) which 
have frequency cut-offs at uj $J 7 2 /& and u> ^ 7/6, respectively. Their contribution to the total radiated energy 
Ed depends on the phase-space in an intricate way, so that the resulting radiation does not have a simple 
universal expression. Specifically, it was found that in the absence of extra dimensions one obtains 

So = C m(^) 2 (^-) 7 3 , Co«8.3xl(r 5 , (5.1) 

with the "basic" relativistic enhancement factor 7 3 . For one extra dimension one has 

Ei =C i m (y) 4 (y) 2 7 3ln 7' d « 1.64 x 10 -8 , (5.2) 
with almost the same (up to the logarithm) enhancement factor. For d 2 one finds 

E d = C d m^f) 7 d + 2 , (5.3) 

with Cd are computed from (|4.10|) . (|4.9|) and (|4. 1 1|) and given above for d = 2, . . . , 6. So the expected enhance- 
ment factor 7 d+2 is regained, with each new dimension adding one power of 7 to the radiation loss. 

Another feature of interest is the spectral-angular distribution of radiation. It was shown that in the usual 
gravity theory without extra dimensions the partial cancelation of local and non-local amplitudes in the case of 
gravitational interaction can be attributed to the fact that in terms of a curved space picture the world lines 
of a massive ultrarelativistic radiating charge and the null geodesic of the emitted radiation stay close to each 
other, so that the formation length of the radiation emitted predominantly in the forward direction is 7 times 
stronger than in flat space. In perturbation theory on a flat background this corresponds to cancelation of 
contributions at high frequencies. In the presence of extra dimensions the emitted radiation is predominantly 
massive from the brane observer point of view, so the trajectories do not stay close together. The resulting 
spectral distribution then has substantial remainder at high frequencies up to 7 2 /6, as illustrated in Figs 4(a)l 

M ■— ' 

The problem considered in this paper is a simplified intermediate step towards the fully gravitational coun- 
terpart, where one is interested in gravitational bremsstrahlung in particle collisions interacting gravitationally 
[24| . The non-local part of the amplitude in this case is due to the three-graviton vertex. The details of this 
problem is the subject of a forthcoming publication. 
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Appendix A: Notation 

1. KK mode decomposition and Fourier transformation - Notation and conventions 

The Fourier decomposition of the bulk fields /imat(cc m , y l ), h' MN (x ,i , y l ), $(x M , y l ) all with periodic conditions 
e.g. h,MN {%> y k + 2ttR) = Hmn{x, y k ) is of the form 

h MN (x,y) = v £ ••• ]T h^ N (x)expU^)= v J2 h MN(xy nkyh/R (AT) 

m — — oo rirf — — 00 ^ ' n 
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Using the representation of the delta-function 

±J2 eink{yk ~ y ' k)/R = sd (y-y')' \ f* inkvh dy d = l[5 nkt0 , (a.2) 

n y k=l 

where V = {2nR) d is the volume of the torus, one obtains the inverse transformation: 

h n MN (*)= I h MN {x,y)e~ m -y k / R d d y. (A.3) 



v 



Four-dimensional fields h r l 1N (x) are then expanded as 

h n MN (x) = 7A4 / e- i{px) h n MN {p) d*P, (A.4) 
where (px)=px = q fl x fl is four-dimensional scalar product, and the final decomposition reads 

hM N (x,y) = tAi^E / hlMP^'^^'^P, (A.5) 

while the inverse transformation is 

h n MN (p) = j " d d y J h MN {x,y)e^- m ^ k l R d i x. (A.6) 

V R* 

Occasionally we will also use another notation for the discrete transversal momenta: p l T — n % / R, i.e. 

h n MN {p) = I h MN (x,y)e^~^yd D x. (A.7) 



with pt y = p T y l ■ From the four-dimensional point of view the zero mode hP MN (x) (n = means all n l = 0) is 
massless, while the n^O modes are massive. Indeed, in the absence of the source term Eq. (|2 . 11 1) reduces to 

d 

(n+ P 2 T )h n MN (x)=o, p! = ^]>> 1 ) 2 ' (a.8) 

i— 1 

where □ = is the four-dimensional D'Alembert operator, while the momenta transverse to the brane give 

rise to the mass term. In the standard scheme [3j one suitably combines polarization modes to get true massive 
gravitons with five spin states for each mass. For our purposes it will be easier to sum over modes using the 
original decomposition. 

When the level spacing is small (e.g. when R^> 6), one can pass from summation over n to integration over 
Pt using 

n v ' 

Here it is implicitly assumed that both the sum and the integral converge. As pointed out in the text, and 
in contrast to quantum Born amplitudes, this is guaranteed in the framework of the classical perturbation 
approach presented here [25| . 

It is worth noting, that upon integration over modes in the case of small level spacing, one obtains the results 
expected in the uncompactified theory in D = 4 + d— dimensional Minkowski space. 

In a similar fashion, expansion of h' MN (x, y) leads to the set of four-dimensional modes h' MN (x), and an ex- 
pansion of the bulk scalar §{x, y) to the set <& n {x). These four-dimensional fields are further Fourier transformed 
to h ,r ^ N {p) and & n (p), respectively. 

2. Integration over angles and frequencies 

In the main text the following integrals over the radiation angle 9 were encountered 

7 sin" 9 

V" = J^—d9, ^ = l-vcos9 (A.10) 
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with integers m,n. 

For 2m > n + 1 one finds to leading order [23 



For n > 2m — 1 one obtains 

2 n-m T (n±l\ T /np, _ \ 

V - = rCn-m + l) 

In the case 2m = n + 1 the leading contribution to the integral is proportional to In 7. For example, one case 
needed in the text was V 2 ~ 4 In 7. 

Calculation of the integrals over the frequency or over the impact parameter involving two Macdonald func- 
tions of the same argument is performed using the formula 36] : 

K,(cz)K v (cz)z^dz = 1 2 j 1 2 1 2 j 1 2 j . (A.13) 

c Q l (a) 



3. Useful kinematical formulae 

The angles in the formulae below are defined in Figure [T] 

u M = 7(1,0, 0,v) , v! = (1,0,0,0) , ip EE 1 — v cos 6 cos ck = 1 — v cos i3 , 

, (ku')b ujb (ku)b ujb , , , 
z = = — , z = = — ip = z yip , 

2y zz' - z 2 - z' 2 = w 2 6 2 sin 2 d , £ 2 = 2 7 zz' - z 2 - z' 2 = uj 2 b 2 sin 2 (9 cos 2 a + b 2 k\ = (yvz' sin df 
— (kb) = £ cos <j) = yz'v sin z9 cos = 72 u cos a sin cos tp = cub sin cos tp 

a , 2 U) 2 b 2 COS $(1 V COS *&) 2 tfi/i i\ fK 1 a\ 

P = yzz' - z 2 = - = 7V -0(1 - VO (A. 14) 



Appendix B: Destructive interference for 7/6 < a; < 7 2 /6 

An alternative proof of the destructive interference effect of the radiation amplitude in the z— region but 
with d < I/7 in higher dimensional Minkowski space will be presented here. In the main text we followed an 
approximation allowing to cover the full angular range. Here destructive interference in the restricted angular 
range will be demonstrated rigorously. 

Start with (|3.20[) . change variable x to ( given by 



£ d( = f(x) dx , f{x) = (z 2 + z' 2 — 272:2:') x + yzz' — z 2 



and integrate by parts twice using 



CkAC) = -K+i(0- 



(B.l) 
(B.2) 



The first integration gives: 



dxe- ix ^K d/2 _ x [<:{x)] 



e -ix(kb) a x —\ 

~^r-K d/2 [C{x)] x=q + I dxK d/2 (()d x 



l —ix(kb) 

1W 



(B.3) 



A second integration by parts leads to 



cr(fc) = \ d 



7 V z 

b d 



yzz' — z^ 



K d / 2 (z)~iq Q 



yzz' — z 2 



L /- . . K d/2+l {z'f 

-\K d/2 {z )-iqi—2 - 

722' V z' — yzz' 



-R 



(B.4) 



where 



q = (kb) - i 



z 2 + z' 2 - 2y zz' 
722' — z z 



v /2 



qi = (kb) - i'- 



2yzz' 



yzz' 



(B.5) 
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and 



R = dx K, 



d/2+l 



(CO*)) 



—ix{kb) 

1W 



fix) 



(B, 



Continuing integration by parts further, one obtains an expansion in terms of go/3 1 and gi(/3 — £ 2 ) . As we 

7 _1 < 1. With this 



discussed before, in the z— region of interest here ip ~ I/7 2 , z ~ 1, z' ~ 7, so that £ 2 ~ /3 
gi ~ 7 and therefore the expansion parameters are: go/3 -1 ~ 7 -1 "C 1, gi(/3 — £ 2 ) -1 ~ 
accuracy one can set go = gi = (kb), ft = jzz' and write: 



z'h 



K d/2 (z')+i^-k d/2+1 (z') 



72^ 



(B.7) 



The first parenthesis in a cancels for v = 1 the leading terms of p (|3.1ip . and the total amplitude j'(fc) = 
p(k) + cr(k) contains only the second parenthesis in (|B.7|) plus the subleading terms mentioned above. Thus, 
the series obtained by integration by parts, converges inside z— cone t? < arcsin7 _1 and establishes the effect 
of destructive interference. 
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